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To accomodate recent RHIC data on 77' multiplicity, we propose a min- 
imal modification of the Witten-Veneziano relation at high temperature. 
This renders a significant drop of rj' mass at high temperature signaling a 
restoration of the U{1)a, and the Goldstone character of 77'. 
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1. Introduction 

On the classical level, QCD with Nj massless flavors enjoys a global 
chiral U{Nf)L U(Nf)ji symmetry. The subgroup with unit determinant, 
namely the SU{Nf)L SU{Nf)ii, is spontaneously broken down to its di- 
agonal part SU{Nf)L^ii = SU{Nf)v- This gives rise to 8 pseudoscalar 
Goldstone bosons: vr's, K's and the r]. 

On the other hand, the non-abelian anomaly in the U{1)a sector, 

= ^^"""FrFr (1) 

explicitely breaks U{1)a symmetry. Coupled with non-trivial topology in 
the field space of QCD, it prevents the spontaneous breaking of U{1)a, thus 
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Fig. 1. The lowest order contribution to the excess of the r?' mass. 

generating the extra mass for 77'. The simplest example of such contribution 
to m^' is the diamond diagram in Fig. [TJ 

Nevertheless, the fate of U{1)a at T > could be drastically changed. New 
RHIC data from central Au+Au collisons [1] on rj' multiplicity shows a drop 
in it's mass by at least 200 MeV inside the fireball 

m^/ (vacuum) = 958 MeV ^ m^/(fireban) = 340l^[}t?^[} ± 45 MeV, (2) 

thus signaling restoration of the Goldstone character of rj' . 

2. Witten-Veneziano relation at finite T 

On the theoretical side, excess in the tj' mass is inferred from the Witten- 
Veneziano relation (WVR) [21 [3] 

m2 + M^, - 2Mi = ^ XYM . (3) 

Jit 

Due to the flavor content of the pseudoscalars, WVR just says that the 
excess in the singlet rj^ mass is coming from the glue sector, i.e. xym is the 
Yang-Mills topological susceptibility. 

Previous studies [4] of WVR indicated, for various lattice forms of xym, 
that 7]' mass increases as xym melts, in contrast to the result ([2]). This 
is happening as T approaches Tch, due to the fact that there f-w{T) starts 
decreasing significantly. 



2.1. Connection between QCD and YM topological susceptibility 

In that light, we propose [5] a minimal modification of WVR first by 
using Leutwyler-Smilga result [6] in the vacuum 
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Fig. 2. The relative-temperature dependences, on T/Tqi^, of x^^^; (99)o^'^i /tt and 
fss- The sohd curve depicts x^^^ for S = 0, and the short-dashed curve is x^^^ for 
6=1. At T = 0, both X 's are equal to xym = (0.1757 GeV)^, the weighted average 
[4] of various lattice results for xym- 



Just like the WVR, Leutwyler-Smilga relation connects the quantities from 
two different theories: x is the QCD topological susceptibility, and Eq. 
shows that it approaches XYM only for large quark masses whereas XYM 
and X are very different for light quarks. S = {qq)o is the condensate in the 
chiral limit. 

We denote by x the whole right hand side of Eq. (HD, and use it at finite 
temperature. 

Importance of this relation comes from the fact that x is driven by the chiral 
quark condensate in the leading order of expansion in small quark masses 



X=-^ + Cm, — =E— • ^ 

Nf ra ^ rrif 



This is the di Vecchia-Veneziano result [6], ^7] . The next term in this expan- 
sion, Cm, is essenital as it keeps xym from blowing up. We fix its value at 
T = by demanding x(0) = Xym- 

2.2. Exploring the thermal dependence 
For the thermal dependence of x we use an ansatz 

s(r)m rs(r) 
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Fig. 3. The relative-temperature dependences, on T/Tc\i, of the pseudoscalar 
masses for 5 — Q. 

This gives 

The interesting window for 5 is then < (5 < 1 since the lower hmit gives 
no thermal dependence for the correction term, and a quadratic one in the 
last equation for x- With 5=1 there is an enhancement of the r]' mass [S], 
and therefore this is the upper limit of interest. 

The interesting window for 5 is then < 5 < 1, since 5 = gives the 
T-independent correction term, while 5 = 1 already leads to precursors of 
the unwanted mass enhancement in the -q' — rj complex. Therefore, 5^1 
is the upper limit of interest, although the results for the T-dependence of 
the meson masses is quantitatively not much different for 6 = 1 than from 
the case 5 = 0, as it is depicted in Figs. [3] and H] respectively. 

3. Results and discussion 

Mesons are constructed as qq bound states via the Bethe-Salpeter equa- 
tion in the ladder approximation. Dynamical quarks are build up from the 
Dyson- Schwinger equation in rainbow approximation. We use the succesfull 
rank-2 separable model [S] for the gluon propagator, which was also used in 
Ref. Rainbow-ladder approximation is the simplest symmetry preserv- 
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Fig. 4. The relative-temperature dependences, on T/Tch, of the pseudoscalar 
masses for S — 1. 



ing truncation, necessary for the correct chiral behavior of the theory. 
The bound state approach enables access only to the non-anomalous part of 
the meson masses, since the ladder Bethe-Salpeter kernel does not include 
diagrams like in Fig. [TJ Therefore, the anomalous part is inferred from Eq. 
([3]). The strategy is to use flavor mass matrices to extract r] and rj' masses 
from the calculated non-anomalous sector. This is presented in detail in 
Ref. m. 

Our main result is presented in Figs. [3] and [H The reduction in the r]' 
mass is around 200 MeV, which is in quantitative agreement with RHIC 
data. This is possible only due to the proposed modification of the WVR 
relation at finite T. Topological susceptibility of pure glue, xyMj is just 
too resistant, with the characteristic melting temperature being Tym = 260 
MeV (see e. g. [21 [TO])- 

In contrast, the pseudocritical temperatures for the chiral and deconfine- 
ment transitions in the full QCD are lower than Tym by some 100 MeV 
or more (e.g., see Ref. [TT]) due to the presence of the quark degrees of 
freedom. In that regard, Eq. ([5]) is essential as it couples chiral restora- 
tion to U{1)a restoration, allowing for x(^) to melt away even sooner than 
f-niT). In the separable model used here we have T^h = 128 MeV which 
is admittedly lower than the so far accepted value around 160-170 MeV. It 
has been shown [T2] that the same model coupled with the gluon degrees of 
freedom in the form of the Polyakov loop cures this discrepancy. 
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